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L INTRODUCTION 

The systems consisting of a two-phase mixture, as clouds of bubbles or drops, have shown many common fea- 
tures in their responses to different external force fields. One of particular interest is the effect of an unsteady pressure 
field applied to these systems, case in which the coupling of the vibrations induced in two neighboring components 
(two drops or two bubbles) may result in an interaction force between them. This behavior was explained by Bjerknes 
[1] by postulating that every body that is moving in an accelerating fluid is subjected to a “kinetic buoyancy” equal 
with the product of the acceleration of the fluid multiplied by the mass of the fluid displaced by the body. The external 
sound wave applied to a system of drops/bubbles triggers secondary sound waves from each component^ of the sys- 
tem. These secondary pressure fields integrated over the surface of the neighboring drop/bubble may result in a force 
additional to the effect of the primary sound wave on each component of the system. In certain conditions, the magni- 
tude of these secondary forces may result in significant changes in the dynamics of each component, thus in the 
behavior of the entire system. In a system containing bubbles, the sound wave radiated by one bubble at the location 
of a neighboring one is dominated by the volume oscillation mode and its effects can be important for a large range of 
frequencies. The interaction forces in a system consisting of drops are much smaller than those consisting of bubbles. 
Therefore, as a first step towards the understanding of the drop-drop interaction subject to external pressure fluctua- 
tions, it is more convenient to study the bubble interactions. 

Effects of the secondary Bjerknes forces are reported in a study concerning cavitation by Kornfeld and 
Suvorov [2], They describe a strange zig-zag motion of the bubbles (“dancing bubbles”) near the surface of an oscil- 
lating piston, driven at high frequencies. The phenomenon received only a qualitative explanation, based on the insta- 
bility of the motion dominated by the inertia forces. Crum [3] studies experimentally the interaction between air 
bubbles immersed in water contained by an oscillating vessel. The vessel is driven at low frequencies (60 Hz) by a 
mechanical shaker. The ambient pressure is reduced below atmospheric conditions in order to increase the natural fre- 
quency of volume oscillations of the bubbles. He provides also an simple model for the relative motion of the bub- 
bles, using a first approximation for the sound field radiated by one bubble at the location of the other one and 
neglecting the inertial forces. His theoretical results compare well with his experimental data at low oscillation fre- 
quencies. However, his algebraic model overpredicts the relative velocity between the bubbles, showing only the 
trends of variation of this quantity. 

The inviscid theoretical approach to secondary Bjerknes forces of Pelekasis [4] and the computational results 
of Pelekasis and Tsamopoulos [5,6] produce useful data for a specific range of frequencies, pressure amplitudes and 
bubble sizes due to the admissible values of Bond and Weber numbers. They predict repulsive secondary Bjerknes 
forces if the frequency of the disturbance is located between the two natural frequencies for volume oscillations of the 
unequal sized bubbles and attraction forces for any other case, including the equal size bubbles interaction. The rela- 
tive motion of the two bubbles is found to be dependent on the frequency of the external acoustic field. When the 

Bond number Bo = (p aR 2 )/a , based on the average acceleration a , the radius R of the larger bubble and the 

properties of the liquid ( p the density, 0 the surface tension coefficient), lies above a critical value, spherical-cap 
shapes are predicted to appear with the deformation confined on the side of the bubbles facing away from the direc- 
tion of acceleration. However, their results are affected quantitatively and qualitatively by the absence of the viscous 
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effects, especially important for the magnitude of the relative velocity and the time in which the two bubbles get to 
collide. This paper presents experimental results and theoretical predictions concerning the interaction and the motion 

of two levitated air bubbles in water in the presence of an acoustic field at high frequencies ( 22-23 kHz ). 

II. EXPERIMENTAL SETUP 

In order to obtain a stable position for the two bubbles within the tank (Figure 1), the dimensions, the quality 
and the amount of water used must be determined in such a way that a stationary sound wave is obtained inside the 
vessel; all these must be correlated with the signal frequency. The exact positions of the bubbles inside the tank are 
determined using an Edmund Scientific laser pointer ML-211 and a milimetric transparent grid. The oscillations are 
induced by a thin wall hollow cylinder piezoceramic transducer made of C-5400 material (Channel Industries) with 

the resonance frequency of 22.5 kHz . The maximum radius of the levitated bubbles with the present amplifier is 

3.5mm . The sinusoidal signal is generated by a Hewlett-Packard 33120A function generator and amplified by a 
Hewlett-Packard 6824A DC Power Supply Amplifier. The signal is checked with a Gould 4050 Oscilloscope (at the 
exit of the amplifier) against any kind of distortions. In these conditions, a stationary sound wave is formed inside the 
tank and the horizontal plane at approximately the middle of water depth corresponds to a pressure minimum. The 
bubbles are levitated at this location with different initial separations. Their motion is recorded and analyzed using a 
Kodak EktaPro 1000 Motion Analyzer, choosing the rate of frames per seconds between 60 and 1000 and adjusting 

the magnification with the Chinon - Hoya Zoom System of lens ( 18 - 108mm ). 

In the principal experiment, two air bubbles of various radii are injected using clean plastic syringes and nee- 
dles at different distances apart and their motion recorded and played back frame-by-frame. Figure 2 gives photos of 
the video screen for a typical experiment, recorded at a rate of 1000 frames per second. The reticle provided by the 
EktaPro system is used to determine the position of the center of mass of each bubble and their radii; a simple calibra- 
tion is performed for each session in order to determine the magnification on the screen. 

In order to properly study their Bjerknes interaction, small bubbles are injected and afterwards their sizes are 
increased gradually until the mutual forces triggered the relative motion. At the high frequency used in experiments, 
it is impossible to measure with accuracy the amplitude of the volume oscillations of the bubbles, this quantity being 
determined indirectly from the equilibrium of forces in the vertical direction. Perfect levitation of the two bubbles is 
obtained easily over the entire central region of the vessel, no significant vertical motion of the two bubbles is 
observed after the initial injection oscillations died down. The quality and stability of the levitation are controlled 

through small changes of the level of the water inside the vessel (around 63 mm ) and the intensity of the acoustic 

field (controlled from the function generator output voltage level with an accuracy of ±1 mV ), the frequency of the 

signal being maintained constant at 22.5 kHz . This paper only presents the attracting bubble results. Other modes of 
bubble interactions will not be discussed here. 

HI. MECHANISM AND DYNAMICS OF BUBBLE - BUBBLE INTERACTION 

Consider now two gas bubbles of nominal radii R 0{ and R 02 , which undergo volume oscillations: 
R { (t) = fl 01 (l + e^incOjl) and R 2 (t) = R 0 2 (1 + e 2 sin(co 2 / + <p)), levitated in an external stationary sound wave 

p A (z, t) = A 0 sin(2rcz/X)sin(coO = Asin(coi) , where A will be a notation for the amplitude of oscillation at a 

fixed location in the liquid column. Volume oscillations of bubble 1 determine the acoustic pressure field around bub- 
ble 2 to be modified with the additional wave (spherical symmetry is assumed) [7]: 

Pi(M) = yi(0^ + 3'i(0[i-(y,(0) 3 K(pt/ 2 .(0)/2) 
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+n - (yi(0)](i/i(0)([pt/ j ,(o - ip A -R x mdp A /(dR x ))vco) 


(i) 


where: y,(r) = (/?,(f))/r , U x (t) = (dR x {t))/{dt) and c 0 =sound velocity in the liquid. The secondary Bjerknes 

force on the bubble 2 due to sound field radiated by bubble 1 will be the average taken over one period of the integral 
of this pressure field on the surface of bubble 2: 

R 12 = (jjPi(r,t)h 3 dS-^ = (j j jVpi(r,t)dV 2 ) (2) 
s 2 v 2 

Assume that the frequencies ct>i and to 2 are both equal with the external applied frequency <n (transient effects 
being very short); also, assume the bubbles having an in-phase oscillation <p = 0 . With these assumptions and by 
neglecting all terms of the order 3 and greater in e y , equations (1 ) and (2) give us the following formulas: 

I'd - * 3£ ;> * ■?]£ ' < 3 > 

where /, j = 1,2 (i * j) . The formulas give a better understanding on the parameters that govern the bubble-bubble 
interaction than the simplified model developed by Crum, whose results showed a pair of equal forces on the bubbles, 
depending directly only with frequency, radii and l/r~ 2 . The validation of our formulas will stand in comparing the 
experimental data (velocities of bubbles) with our model predictions and with Crum’s experiments and model results. 

The forces which are considered in studying the motion of the bubbles in the horizontal plan are the secondary 
Bjerknes forces given by (3) and the drag forces: F Di \ = (pv I ^5^C £)( )/2, where Re ( = (pv l (2i? 0l ))/2 is Reynolds 

number for bubble i based on the instantaneous velocities v- and nominal radii R 0i , and Sj = is the frontal 
area of the bubble. The drag coefficient is determined by the model proposed by Moore [8]: 
C Di = 48/te r ~ 1 (l -2.2 Re ~°‘ 5 ) . The forces generated by the gradient of pressure amplitude in the horizontal plane 

are neglected as the hydrophone measurements indicated a smooth distribution of the amplitude A .The differential 
equations governing the motion of the two bubbles are nonlinear, of second order with respect to time: 

m t r i = |F^| - \F Di \ , where m l is the induced mass of bubble i which is equal to the half of the mass of the dis- 
placed liquid).This relation is transformed by changing the independent variable to r(t) = r { ( t ) - r 2 (t) and choos- 
ing the unknown functions to be the velocities of each bubble, v ( . The amplitude of the pressure signal will be 

considered the same around each bubble: A { = A 2 = pco 2 R^t\ = p 0) 2 /?^£| since co$ « co 2 for the investigated 
frequencies. Further substitution results in the following system: 

(v 1 + v 2 )v/ = G t r~ 2 - ((72u)/ /?&) v f [ 1 - 2.2(v/ (2v,/? oi )) 0 ' 5 ] 

where: G; = 1.5R 0 3 ; (o 2 [e y (e ; + 2e i )] , i,j = l,2(i*j), v{ = (dv x )/(dr), and v 2 ' = ~(dv 2 )/(dr) . The 

length scale and the velocity scale used to non-dimensionalize (4) are: L = R x and U = (a/(p/?,)) 0 5 , where the 
larger bubble (index 1) is considered to be the left one. This nonlinear system is solved numerically with the initial 
conditions imposed by the initial velocities measured during each experiment: v x (r=d 0 ) = v 10 and 

v 2 (r=r/ 0 ). = v 20 . The solution of the system is used to compute and plot the relative velocity of the bubbles 
v(r) = v,(r) + v 2 (r). 
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IV. RESULTS AND DISCUSSIONS 


Figure 3 presents the results obtained by applying the model proposed by Crum and our model to our experi- 
mental data in order to compare the accuracy of both approaches. The comparison shows clearly that the simplified 
dynamics considered by Crum overpredicts the relative velocity even at large bubble separation. Neglecting the iner- 
tia forces and using a first order expression for the acoustic field radiated by one bubble in determining the secondary 
Bjerknes forces formulas are the principal causes for its limited accuracy. On the other hand, the new analytic model 
presently proposed gives results in good agreement with the experiments until the dimensionless distance between 

bubbles decreases below a minimum value: r m = 3 as shown in Fig. 3. The differences between theory and experi- 
ment become significant for r <r m , where the two neglected effects mentioned before (drag modification and acous- 
tic interference) are important. At the moment of collision, secondary Bjerknes forces attain the maximum 
magnitude, which is less than 40% of the buoyancy on the larger bubble. Based on the proposed model, we designate 
two new dimensionless parameters which determine the time in which the bubbles reach the impact 

point: T = d$/(R 0l R 02 ) (containing the geometrical conditions of the interaction) and a = A/(P 0 (z)) (related to 
the acoustical conditions). Figure 4 presents the experimental data translated in dimensionless variables. These plots 
of the coordinates of center of mass of each bubble in time show that: At ~ P’or 1 with positive values for n and s . 
The average values correlated from (he experimental data are: n = 4/3 and r « 1/2 . 
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Distilled water 





New Differential Model 



Fig. 3. Relative velocity of the bubbles. Fig. 4. History of x-coordinates of the centers 

Comparison of new experimental data, of mass of each bubble a between 0.30 

Crum's model and new differential and 0.49 

model: T=390.6 and a=0.49 




